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Abstract. The overdamped motion of a Brownian particle in randomly switching piece-wise metastable 
linear potential shows noise enhanced stability (NES): the noise stabilizes the metastable system and the 
system remains in this state for a longer time than in the absence of white noise. The mean first passage 
time (MFPT) has a maximum at a finite value of white noise intensity. The analytical expression of MFPT 
in terms of the white noise intensity, the parameters of the potential barrier, and of the dichotomous noise 
is derived. The conditions for the NES phenomenon and the parameter region where the effect can be 
observed are obtained. The mean first passage time behaviours as a function of the mean flipping rate of 
the potential for unstable and metastable initial configurations are also analyzed. We observe the resonant 
activation phenomenon for initial metastable configuration of the potential profile. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 02.50.-r Prob- 
ability theory, stochastic processes, and statistics - 05.10.Gg Stochastic analysis methods (Fokker-Planck, 
Langevin, etc.) 



1 Introduction 

The thermally activated escape from a metastable state in 
fluctuating potential is of great importance to many nat- 
ural systems, ranging from physical and chemical systems 
to biological complex systems. A particular challenging 
direction are systems far from thermal equilibrium, ei- 
ther due to non-thermal noise or external deterministic 
periodical forces [1121314151617181911011 lj . A typical prob- 
lem is the enhancement of stability of metastable and un- 
stable states due to the external noise |7l8l9ll(JfrT| . The 
noise-enhanced stability (NES) phenomenon was observed 
experimentally and numerically in various physical sys- 
tems (see Refs. |5I6I7I8I9I10I11I12I13I14| and, as a recent 
, review, Ref. HS1). The investigated systems was subjected 
to the action of two forces: additive white noise and regu- 
lar force. The regular force was fixed or periodical in time, 
so that metastable state appeared for a half of period and 
unstable state appeared for an other half. The noise en- 
hanced stability effect implies that, under the action of 
additive noise, a system remains in the metastable state 
for a longer time then in the deterministic case, and the 
escape time has a maximum as a function of noise inten- 
sity. We can lengthen or shorten the mean lifetime of the 
metastable state of our physical system, by acting on the 
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white noise intensity. The noise-induced stabilization of 
one-dimensional maps |16I17| . the noise-induced stability 
in fluctuating bistable potentials ^H], the noise induced 
slowing down in a periodical potential 1 9 20121] , the noise 
induced order in one-dimensional map of the Belousov- 
Zhabotinsky reaction |22I23| . and the transient proper- 
ties of a bistable kinetic system driven by two correlated 
noises [21], are akin to the NES phenomenon. Moreover, 
this effect is at the basis of resonant trapping 125]. Even 
though the previous theoretical papers analyzed NES phe- 
nomenon in systems with fixed or periodically driving 
metastable and unstable states, the model of randomly 
switching metastable state is more realistic in many cases, 
e. g. when we describe the generation process of the carrier 
traps in semiconductors. Despite its experimental impor- 
tance, the theory of fluctuating barrier crossing is not well 
developed for arbitrary noise intensity. In the present pa- 
per we obtain and study analytically the NES effect in a 
system described by a potential, which randomly switches 
between metastable and unstable configurations. We de- 
fine the lifetime of metastable state as the mean first pas- 
sage time (MFPT) and obtain the conditions when it can 
grow with white noise intensity. The mean first passage 
time behaviour as a function of the mean flipping rate of 
the potential for unstable and metastable initial config- 
urations is also analyzed. We observe the resonant acti- 
vation phenomenon for initial metastable configuration of 
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the potential profile, and monotonic behaviour for initial 
unstable configuration. We can describe therefore, with 
the same theoretical approach, two noise-induced phenom- 
ena |7l1 ttM<MiX\ . 



2 The model 

We consider one-dimensional overdamped Brownian mo- 
tion in a randomly switching potential profile U (x) 



dx 
~db 



= f(x)-ar,(t) + at): 



(1) 



where £ (t) is the white Gaussian noise with zero mean and 
+ T )) = 2Z?<5(r), 77(f) is a Markovian dichotomous 
process, which takes the values ±1 with mean flipping rate 
v, and f(x) = —dll(x)/dx. We investigate the mean first 
passage time with the reflecting boundary at the point 
x — and the absorbing boundary at the point x = b 
(b > 0). The calculation technique for the MFPT of non- 
Markovian process x(t) without white noise (£(f) = 0) has 
been originally developed in Ref. [2H] and, then, general- 
ized by various authors (see, for example Refs. 



The exact equations for mean first passage times, which 
take into account both the dichotomous and the white 
noise terms in Eq. (JTJ, were derived in Ref. [33], where 
the authors solve the delicate problem to construct correct 
conditions at the reflecting boundary, since all previous 
works dealt with absorbing boundaries only. Nevertheless 
the cited boundary conditions (3.7a) and (3.7b) in Ref. 
|33| are only valid for the special unlikely case of immedi- 
ate reflection. For our purposes we use the same conditions 
at the reflecting boundary as in Ref. [32]. Thus from the 
backward Fokker-Planck equation we obtain the coupled 
differential equations for the MFPTs in our system 
(see Ref. [33] and Appendix A) 



DT'l + [/ (x) -a]T' + + v (T_ - T+) = -1, 
DT'l + [/ (x) +a]T'_+v (T+ - T_) = -1. 



(2) 



Here T+(x) and T—(x) are respectively the mean first 
passage times of the boundary x — b for initial values 
77(0) = +1 and 77(0) = — 1, with the Brownian particle 
starting from the initial point x (0 < x < b). Our condi- 
tions at the reflecting boundary x = and the absorbing 
boundary x = b are (see Appendix B) 



T4(o) = o, r±(6) = o. 



(3) 



Let us introduce for convenience two new auxiliary func- 
tions 



T = 



9 



T_ -T, 



(4) 



2 2 
and rewrite Eqs. in more simple form 

DT" + f (x) T' + a©' = -1, 
DO" + f (x) 0' + aT' - 2v0 = 0. (5) 

The boundary conditions for the functions T(x) and O(x) 
follow from Eqs. © and Q 

T' (0) = 0' (0) = 0, T (b) = O (b) = 0. (6) 




Fig. 1. Two configurations of fluctuating potential U(x), cor- 
responding to the metastable and unstable states. 

Equations (J5J with the boundary conditions © were solved 
in Ref. [33] for a model with constant force 



fix) = const. 



< x < b. 



(7) 



In the present paper we consider more general model with 
step- wise force (or piece- wise potential U{x)) as 



fix) 



dU(x) 
dx 



k ■ 9(x - L), 



< L < 6, 



(8) 



where 9(x) is the step function and is a constant. As we 
can see from Fig. [fl if < a < k we have a metastable 
state for 77(f) = +1 and an unstable state for 77(f) = — 1. 
In particular cases when L — > or L — > b the force 
coincides with JJJ. After removing T{x) from Eq. (J5J 
we obtain third-order linear differential equation for the 
variable 0(x) 



where 



2/ (a) 
D 



0" + F(x)0' 



M (x) 

D 2 



F[x) 



P(x) , / jx) 
D 2 



D 



■ 7 



7 



D-^D- (10) 



We will consider the mean first passage times T±{0) with 
the starting position of Brownian particle at x = 0. In 
the absence of switchings and white noise the dynamical 
escape time T + (0) for the initial metastable state is equal 
to +00 and the dynamical escape time T_ (0) for the initial 
unstable state equals 



b-L 



(11) 



We solve Eqs. (JEJ and J3J for regions < x < L and 
L < x < b separately. First of all we find the solutions in 
the interval (L, b) 



i=i 
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T (x) = C4 



X — L 



[ e M*-L) - 1] 

k + X,D 



(12) 



where Cj (i = 1 -j- 4) are unknown constants and Ai, A2, 
A3 are the roots of the following cubic equation 



A (AD + fc) - T 2 A - 2vk = 0, 



(13) 



where r = 7D. Using graphical representation it can be 
easily shown that algebraic Eq. (|13fl has three real roots: 
one positive and two negative ones. 

Taking into account the conditions I© at the reflecting 
boundary i = 0we obtain for region (0, L) 

O (x) = C5 ( 1 + — jj- cosh 72? J + (D sinh 7 a; — fx) , 



T(x) 



(cosh 7a; — 1) 



p2 p4 

C5 (D sinh ■yx — rx) + eg . 

ra 



From Eqs. 10} and lfH)l we get finally 

r±(o) =t(o)t0(o) = c 6 t4 



(14) 



(15) 



The six unknown constants Cj (2 = 1 -j- 6) can be deter- 
mined from the conditions JSJ) at the absorbing boundary 
x = b and the continuity conditions of the functions 0(x), 
0'(x), T(x), T'(x) at the point x = L. Thus, we obtain 
the following system of algebraic linear equations for the 
constants & 



3 

£ 



3 

£ 



c t a (e"' - 1) 
fc + XiD 



b-L 



C4 



2*aD 

c, = cn [ 1 H 5- cosh 7L 



(£> sinh 7 L-rL) + ^, 



^ c,:Ai 



i=l 
3 



yy '7 A," 



2i/r a , 

C5 — 5- sinh 7 L + — (cosh 7 .L — 1) . 
1 z 

2vL a 2 



k + \ x D r 2 r 3 
1 

a A: ' 



sinh 7L 



i=l 

H c 5 (cosh7L — 1 



Ci = c 6 - -j^ Yi (cosh 7 L - 1) 

-c 5 ^(Dsmh~/L-rL), 
1 a 



(16) 



where \ii — Xi(b — L). The equation i|15|) with algebraic 
system (|16fl is the exact solution of Eqs. for the force 
||S}. In the limits L — ► or L — ► b, the expression i|15f) 
coincides with that obtained in Ref. [2] for the force Q. 



3 NES phenomenon conditions 

Our aim is to investigate the noise enhanced stability ef- 
fect. In other words we are interested in the situation, 
when the escape time grows with noise. We expect to find 
this growing in the limit of small intensity of the white 
noise. Hereafter we write the mean first passage times 
T±(0) as T±(y,D) to point out the dependence on the 
parameters v and D of two noise sources. On cumber- 
some rearrangements (see Appendix C), from Eqs. (|15|l 
and l|16|l we obtain in the limit D — > 

T-(u, D) ~ T_ (u, 0) + ^G(q, w,s) +o(D) , (17) 



where 
T_(i/,0) 



2L vL 2 b 



L q(l-q) 



(1 



a a* k 2v 
is the MFPT in the absence of white noise, and 

„3^ 



(18) 



G{q,u,s) 



q 2 ) se 



5<T - 5q 3 



(l + q)(l-q 2 ) 2(l + q)(l-q 2 ) 
q (1 + q - 5q 2 - 3q 3 - 2q 4 ^ 



2(l + g ) 



1 



2lo (3q 2 



q(l-q 2 ) 



2to 2 
Q 2 ' 



Here q, cj and s are the dimensionless parameters 



a vL 2lj 

^ k ' fe ' 1 — q 2 



b 

L ~ 1 



(19) 



(20) 



As one would expect, the MFPT l|18(l coincides with the 
nonlinear relaxation time (NLRT) for the same system 
(see Eq. (32) in Ref. |3S]) because in the absence of dif- 
fusion (D = 0), after crossing the point x — b Brownian 
particles cannot come back in the interval (0, b) again. 

The sign of the term G(q,uj,s) in Eq. (|19fl allows to 
obtain the conditions to observe the NES effect in the 
system investigated. We can write these conditions as 



G(q,u,s) > 0. 



(21) 



Let us analyse the structure of NES phenomenon region 
on the plane (q,oj). In the case of very slow switchings 
v — > (lu — > 0, s — > 0) Eq. (J^T} takes the form 



us < 



q(l-q) (5q 3 + 5q 2 -q-5) 



2[6-2q-q 2 (b/L + 11) + 2q 3 + 3g 4 (b/L + 1)] ' 

q> 0,8024. (22) 
In the case of q ~ 1 we obtain from Eq. (|21|) 

a (1 — q) 15, . 1 
W< T^TT' 2 + 2( 1 -^ < - < 2TT^)' (23) 

where xo is the positive root of algebraic equation 

e x = x + 2 
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Fig. 2. The shaded areas are the regions of the plane (q, ui) 
where the NES effect takes place. The parameter is b/L = 1.5. 

approximately equals xq — 1.1463. The two shaded areas 
on the plane (q,w), where the NES phenomenon takes 
place, are shown in Fig. Both noise enhanced stability 
areas connect at the parameter b/L < 1.2655. 

Now we may conclude that the NES effect occurs at 
the values of q near 1, i.e. at very small steepness k — a = 
k(l — q) of the reverse potential barrier for the metastable 
state (see Fig.^l. Only in such a situation, Brownian par- 
ticles can move back into potential well in the interval 
(0, L) from the region L < x < b, with a low noise inten- 
sity. In Ref. JOj w e obtained the parameter region of NES 
effect for periodical driven metastable state. This region is 
larger than that obtained here for randomly switching po- 
tential. In our opinion, this narrowing of NES phenomenon 
area is attributable to the uncertainty of the first switch- 
ing time. Moreover, because of the exponential probability 
distribution of the time interval between jumps of the di- 
chotomous noise, which randomly switches our metastable 
state, the average escape time has a large variance and 
thus the analogy between the two cases falls. 

It must be emphasized that noise enhanced stability 
regions of the mean first passage time are inside that ob- 
tained for the nonlinear relaxation time with the same pa- 
rameter b/L (compare Eq. I|19|) with Eq. (31) in Ref. |36|) 
because the nonlinear relaxation time takes into account 
a repeated reentries of Brownian particles in the inter- 
val (0,6). The NES effect disappears when the absorbing 
boundary is placed at the point x = L. In fact by putting 
b — » L, k — » +oo (q — > 0, s — > 0) in Eq. (|17|) we obtain 

, . 2L vL 2 D (h &vL 2v 2 L 2 \ , s 

T_0 ~— + — ^ - + + —a- ■ 24 

a a A a A \2 a a 1 ) 

The behaviours of the MFPTs T_(i/, D) for the initial 
unstable configuration, normalized to the value obtained 
in the absence of white noise T_(^, 0), vs the white noise 
intensity are shown in Figs. E3 and 0] for the two NES areas 



T-(v,D) 
7.(v,0) 




Fig. 3. The normalized MFPT vs the white noise intensity 
D, for three values of the dimensionless mean flipping rate 
uj = (uL)/k: 0.01 (curve 1), 0.04 (curve 2), 0.07 (curve 3) 
(lower area in Fig. |3J. Parameters are a — 0.9, k = 1, L = 1, 
and 6 = 1.5. 



1 .0003 



7"-(v.P) 
T_(v,0) 




Fig. 4. The normalized MFPT vs the white noise intensity 
D, for three values of the dimensionless mean flipping rate 
uj = (uL)/k: 0.7 (curve 1), 0.8 (curve 2), 0.9 (curve 3) (upper 
area in Fig. Parameters are a = 0.97, k = 1, L = 1, and 
b = 1.5. 

represented in Fig. |21 As we can see the effect decreases, 
for fixed q, when u> decreases in the upper area of Fig. [3 
and when ui increases in the lower area of Fig. [3 At the 
same time, the NES phenomenon for the upper area in 
Fig. [2 is very small (see Fig.0J. 

Now we analyse the behaviour of the escape times 
T-{y,D) and T + (u,D) as a function of the mean flip- 
ping rate v. From Eqs. I|15|l and l|16|) we obtain these 
behaviours, which are shown in Fig. f° r a fixed value 
of the white noise intensity D — 0.1. 

We find the resonant activation phenomenon :26 28 
for T + (y,D) and monotonic behaviour for T-{y,D), In 
the same Fig. [3]we report the arithmetic average of these 
quantities: f(v,D) = [T + (v,D) + T-(v,D)]/2, which is 
the average escape time analyzed by Doering and Gadoua 
in Ref. Starting from the metastable configuration, 
the mean first passage time T+(v, D) decreases with the 
mean flipping rate, because the potential fluctuations in- 
duce crossing events of the potential barrier. On the other 
hand, starting from the unstable configuration, the av- 
erage escape time T-(v,D) increases monotonically be- 
cause the potential fluctuations stabilise the initial unsta- 
ble state. For very fast potential fluctuations, Brownian 
particles " forget" the initial configuration of the potential 
and all the mean first passage times (7+(z/, D), T_(y,D), 
and T(y, D)) join in the same asymptotic value: the av- 
erage escape time required to cross the average barrier 

EE3 
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Fig. 5. Logarithmic plot of the mean first passage times: 
T+(v,D) (initial metastable configuration), T-(v,D) (initial 
unstable configuration), T(v,D) (arithmetic average), as a 
function the mean flipping rate v. Here D = O.f, a — 0.5 
and the other parameters are the same as in Figs. [3] and [I] 



4 Conclusions 



We have investigated the noise enhanced stability (NES) 
phenomenon in one-dimensional system with Gaussian ad- 
ditive white noise and a potential randomly switched by 
a Markovian dichotomous process. We have derived gen- 
eral equations to calculate the mean first passage times 
with one reflecting boundary and one absorbing bound- 
ary. For piece-wise linear potential we obtain the exact 
solution for MFPTs as a function of arbitrary white noise 
intensity and arbitrary mean flipping rate of the potential 
barrier. By analysing the derived equations we obtain an- 
alytically the region of the NES phenomenon occurrence. 
We find that this effect can be observed only at very flat- 
tened sink beyond the potential barrier, i.e. at the real 
absence of the reverse potential barrier in the metastable 
state. Moreover we note that the noise enhanced stability 
phenomenon is related to rare events, that is when some 
particles are trapped in the metastable state. Therefore 
to better understand this effect it should be interesting 
to analyse the probability distribution of the escape time. 
This will be done in a forthcoming paper. By investigat- 
ing the behaviour of the mean first passage times as a 
function of the mean flipping rate of the potential we find 
the resonant activation phenomenon for initial metastable 
configuration and monotonic behaviour for initial unstable 
configuration. With our theoretical approach we are able 
to describe, therefore, two different noise-induced effects. 
Our exact analytical results, concerning the enhancement 
of stability of metastable states for one-dimensional sys- 
tems, can be a good starting point to extend and improve 
the same theoretical apparatus to more complex systems 
with fluctuating potential barriers. 
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5 Appendix A 

Let us consider a stochastic Langevin equation with two 
random forces 



x = f(x)+g (a:) r,{t) + h (x) £ (t) 



(25) 



where £ (t) is the white Gaussian noise with zero mean 
and correlation function (£ (t) £ (t + r)) = 2D 5 (r), r\ (t) is 
a Markovian dichotomous noise with flippings mean rate 
v and values ±1, / (x), g (x), and h (x) are arbitrary func- 
tions. To derive the closed equation for the joint probabil- 
ity density of the random processes x(t) and rj(t) we start 
from the following expression 

W [x, y, t) = (S(x- x(t)) S (y - r,(t))) (26) 

and then apply well-known method based on functional 
correlation formulae (see, for example, Refs. |35I36| L 
Upon differentiation of Eq. I(26|) on i, we obtain 

dW d 

(x(t)6(x-x(t))S(y-rj(t))) 



dt 



dx 

+ (S(x-x(t))^S(y-r,(t)) 



(27) 



Substituting x (t) from Eq. (|25|l and taking into account 
the definition 1)26(1. we can rewrite Eq. 1(27(1 as 



dW 
~dt 
d_ 

dx 



d_ 

dx 



[f(x)+yg(x)]W 



h(x) (Z(t)6(x-x(t))S(y- V (t))) 







+ (S(x-x(t))-S(y-r 1 (t)) 



(28) 



To split the functional average in Eq. (|28|l we use Furutsu- 
Novikov's formula for the white Gaussian noise £ (t) 37 38 



(£(i)F t [£])=£> 



SFt j£j 
8t (t) 



(29) 



where F t [£] is an arbitrary functional depending on the 
history of random process £ (t) (0 < r < t). Replac- 
ing F t [£] with 5{x- x (t)) S (y - rj (t)) in Eq. ^ and 
taking into account that, in accordance with Eq. 1)25(1 . 
Sx (t) /S^ (t) — h(x(t)), we obtain 

(£ (t) S(x-x (t)) 8(3,-1 (t))) = -D^-h (x) W. (30) 

Further let us rearrange the last part of second average 
in Eq. I(28() . Using the translation operator and taking into 
account that for dichotomous noise with the values ±1: 
v 2k (t) = 1, i] 2k+1 (t) = r] (t) (k = 1, 2, . . .), we arrive at 

p (y -*{t)) = | exp 

-|sinh L(t)-^)S(V) = -77(t)sinh ^ ) ,>(.,/). CU ; 



Mt) Ty )S(y) 
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Substituting Eqs. (jSOJ) and in Eq. we obtain 

dW d d d 

~df = ~~dx + V9( x )] W + D—h(x) — h (x) W 



f)(t)S(x-x(t))shxh [ — ) S(y) 



(32) 



Now we use the Shapiro-Loginov's formula for a Marko- 
vian dichotomous noise r\ (t) |39| , to express the functional 
average in Eq. (|32[1 in terms of the joint probability dis- 
tribution 



-2v (v(t)Rt[r))) 



(33) 



where Rt [rj\ is an arbitrary functional depending on the 
history of random process r] (t) (0 < r < t). As a result, 
for the functional average in Eq. (|32|) we have 



v(t)6(o 



x(t))sinh[ — )8(y) 



-2v U(x - x{t)) sinh I ry (t) - 5(y) 

-v {S(x - x(t)) [S(y + r)(t)) - S(y - V (t))}) . (34) 

Substituting Eq. (|34(l in Eq. (|32J) and taking into ac- 
count Eq. (|26|l we obtain finally the following forward Kol- 
mogorov's equation for the joint probability distribution 



[/ (x) + yg (x)} W + D 



dW _ d 
dt dx 

+ v[W(x,-y,t)-W(x,y,t)} 



d_ 

dx 



h{x) 



W 



(35) 



According to Eq. I|35l) we can reconstruct the back- 
ward Kolmogorov's equation for the conditional proba- 
bility distribution W(x, y, t \xo,yo,to) of the Markovian 
vector-process {x (t) ,rj(t)} 



dW dW 

"W = Lf ( x o) + Vo9 (xo)] jr- 
oto axo 



D 



h {x Q ) 

OXQ 



W 



+v [W(x, y, t \x , -y , t ) - W[x, y, t \x , y , t )] . (36) 

Since the conditional probability density depends only on 
the difference r = t — to , Eq. (|3*6*|l can be rewritten as 



dW dW 
-g^r = [f (xo) + yog {x )] 



D 



h(x ) 

dx 



W 



+v[W(x,y,T\x ,-y ,Q) - W(x,y,T\x o ,y o ,0)] . (37) 

Let us introduce the random first passage time 9. The 
probability Pr{# > r} = P (xo,yo, T ) that random pro- 
cess x (t) remains between the absorbing boundaries x = 
L\ and x = L2 during the time interval (0, r) is 



P{xo,Vo,t) 



dx 

Li — \ 



W(x,y,T\x ,y ,0)dy. (38) 



Taking into account Eq. I|38[) and integrating Eq. (|37|) on 
x and y we find 



dP dP 
+v [P(x , -yo, t) - P(x , y , t)] . 



h(x ) 



d 
dxn 



P 

(39) 



To obtain the probability density of first passage time we 
should differentiate the probability Pr{6> > t} on t 



w(x a ,y ,T) = 



dPi{d<r}_ dP{x ,yo,T) 



dr dr 
As a result, we have the same Eq. (|3^)l for w (x Q , y , r) 



Ow dw 

-s- = [f ( x o) + yog (x )} -t^ + d 

OT OXq 

+v[w(x ,-yo,r) - w(x a , yo ,T)} . 
For the mean first passage time 



h (x ) 

OXq 



$ (xo,yo) 



Tw(x ,y ,T)dT 



(40) 



(41) 



we obtain from Eq. 14U|) the following differential equa- 



tion 



D 



h(x ) 



d 



dx 

+v [d{xQ, -y ) - ti(x , y )] 



$ + [f (xo) + y g (x )} 



dx 



(42) 



Since yo = ±1, Eq. (|42(1 is equivalent to the following set 
of equations for times T + (xo) = ${xo, +1) and T_ (xq) = 
i?(xo,-l) 

dT 4 



L X0 T + + [f (x ) + g (xo)} -j^ + v(T-- T+) = -1, 

u(T + -T_) = -l, (43) 



dxo 



L X0 T_ + [f (x ) ~ g (xq)} 

dxo 

where we introduced the operator 



L n 



D 



h(xo)- — 
dxo 



(44) 



Equations H43|l are similar to well-known equations for 
mean first passage times obtained in Ref . 1331 . But authors 
used the fto's interpretation of Langevin equation (|25|l 
and, as a result, obtained different equations. For h (x) = 
f, Eqs. H43fl coincide with the equations derived in Ref. 
PSj. Finally, if we put in Eqs. (jUII and (jUJ): g (x) = -a, 
h(x) = 1, we obtain Eqs. @. 



6 Appendix B 

We explain the origin of conditions © for the mean first 
passage times at a reflecting boundary. For sake of def- 
initeness, we consider a left reflecting boundary at the 
point A. Let Brownian particle reaches the reflecting bound- 
ary at the instant £o — Ato- For the next small time period 
Ato, dichotomous noise r](t) can switch to the opposite 
state with the probability P(Ato) or remains in the initial 
state with the probability (f — (3(Ato))- The scenario of 
Brownian particle behaviour for this time interval is the 
following: it transfers to the point A+Axo (Axo > 0) with 
some fixed probability a or stays at reflecting boundary 
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with the probability (1 — a). We denote the probability a 
as ev_ for the case when dichotomous noise switches and as 
a + for the opposite case. As a result, the joint conditional 
probability distribution reads 

W(x , z , t \A, y , t -At ) = (3{At )S (z + y ) ■ 

[a-5 {xq-A- Axo) + (1 - a-) S (x - A)} 

+ [1 - P(At )} 5 (z - y ) [a+S (x -A- Ax ) 

+ (1 - a + ) 5 (x - A)}. (45) 

Substituting Eq. (|45|l in Smoluchowski equation 

r + co r- + oo 

W(x,y,t\A,y ,t - At ) = dx dz 

J A J — oo 

W[x ,Zo,t \A,y ,t - Ato)W(x,y,t\xo,zo,to) 
we obtain 

W(x,y,t\A,y ,t - At ) 

= a-f3{At )W(x, y,t\A + Ax , -y , t ) 

+a+ [1 - P(At )} W(x, y,t\A + Ax , Vo, to ) 

+(3{At ) (1 - a-) W{x, y, t \A, -y , t ) 

+ (1 - a+) [1 - [3(Ato)} W(x, y, t \A, y , to). (46) 

By expanding the conditional probability distributions in 
Eq. I(46j) in power series in Ato, up to a linear terms, and 
in Axq, up to a quadratic terms, we find approximately 



-At 



dW+ 



dto 

+a-/3{Ato)Ax 



P{At ) (W- - W+) 



( 



dW- 



V dxo 
+a+Ax [1 - P{Ato)] 



Axo d 2 W- 



A 

dW+ 



dxo 



dxl 



Axo d 2 W+ 



dxl 



(47) 



where for simplicity we introduce the notations W+ = 
W{x,y,t\A,y ,t ) and W- = W(x,y,t\A,-y ,t )- 

For Markovian dichotomous noise we have f3(Ato) = 
vAto- After substituting this value in Eq. 1|47[) and divid- 
ing by u+Axo both sides of the equation, we arrive at 



At dW+ 



vAto 



a+Axo dto a+Axo 



+ (1 - vAto) 



dxn 



{W- - W+) 
Axo d 2 W + 



dxl 



+vAt 



a_ 



dxn 



Axo d 2 W. 



dxl 



(48) 



Since for a diffusion process x(t): Axo • 
Eq. we obtain in the limit At — > 



dW + 



dx 



s/DAt , from 



(49) 



By differentiating Eq. (|41|l (see Appendix A) and using 
Eq. gUJ) we find 



T[ (A) = 0, T'_ (A) = 0, 



(50) 



that are the boundary conditions used in our paper. 

If we request an immediate switching of dichotomous 
noise rj(t) when Brownian particle reaches the reflecting 
boundary A, we must put (3(Ato) = 1 in Eq. 147(1 . As a 
result, we have 



-At 



dW+ 
~dlo~ 



a-Axo 
W. 



( dW_ 



Axo d 2 W_ 



V dxo 
W+. 

In the limit Ato — > we obtain from Eq. I(51|l 

W+ = 



dxl 



(51) 



(52) 

Taking into account equality (|52() in Eq. i|51[l and dividing 
both sides of this equation by a~Axo, we find in the limit 
Ato 

dW- 



dxo 



0. 



(53) 



Equations (|52() and (|53|l are equivalent to the following 
conditions for MFPTs at the reflecting boundary 



T + (A)=T_(A), T'(A) = 



(54) 



and 



T+ (A) = T- (A) , T'_ (A) = 0. (55) 

It would be emphasized that conditions l|54|l and l|55|l were 
previously derived in Ref. 33 by more complex procedure. 
We must choose the conditions (|54|l or, alternatively, the 
conditions l|55|) from physical considerations. 



7 Appendix C 

Using Eqs. (pU we can eliminate the unknown constant 
Cq from Eq. (|15[1 and after rearrangements obtain 



b vL 2 
k + T 2 



t ±(0) = t: + — + E ac * 



2=1 

r 2 2 V l r 2 
.___ C5 _ TC5 



k + \ t D + T 2 



i f a\ 



(56) 



The unknown constants Ci , 02,03, C5 involved in Eq. (|56|l 
must be found from the following set of equations (see 
Eqs. J32J) 



2vV 



E c * e ' 

i=l 

A / 2vD 

c i - c 5 ( 1 H — cosh 7L 

i=l ^ 

<xD . a aL 
= — t sinn 7L H , 

r 3 ' 2vk r 2 ' 

QiisT^ (x 
} J CiXi - c 5 — — sinh 7L = — (cosl^L - 1) , 



i=l 
3 

£ 



c 5 ^ (cosh7_L - 1) 



2uL a , 1 
_ + _ slnh7L __. 



(57) 
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To investigate the NES phenomenon, we look for asymp- 
totic expressions of Eqs. (|56|l and (|57|l for D — > 0, with an 
accuracy of linear terms. To do this we need approximate 
expressions of the cubic equation roots. From Eq. (|13|) we 
have 



Ai + A2 + A3 — 
2vk 



2k 



Ai A2A3 



(58) 



and, as a consequence, we can put the roots in the follow- 
ing form 



A 2 + A 3 D, 



(59) 



with the unknown constants A±, A 2 , A 3 . Substitution of 
Eq. (Si in Eq. JT3J) gives 



Ai 
A 2 

A3 



v 



D 

k — a 

D 
2vk 



k + a 
v 



k — a 

2vD (fc 2 + a 2 ) 
(fc 2 - a 2 ) 2 



(60) 



(7 



As it is seen from Eq. IjtjUfl . in the limit D 
+00) we have: fi± — > —00 and \i% — > —00. Therefore, the 
negligibly small terms with e Ml and e M2 can be neglected 
in Eq. I|56l) and in the first equation (|57|l . and we find 



C3 



2vk 



(61) 



Using the following approximate expressions 

P 7i 



sinh 7L ~ COSI17L 
we find C5 from Eqs. (|57|l 



(7 



-co 



C5 



2vr 3 ' 



(62) 



Substituting Eqs. JHTJ and (fr!2*)l in Eq. f5H|l . and solving 
the system l|57l) . we obtain 



T±(0) 



r 2 



6 

fc ' r 2 2vk{k + \ 3 D) 
a 2 [k{r 2 - AvD + 2vLF) - r 2 {F + DX 3 )} 

r i x 3 {k + DX 3 )(k + D\ 3 + r) 

a 2 e~^ [(fc + DX 3 )(k + 3DA 3 ) - F 2 ] 

2vkF{k + DX 3 ){k + DX 3 + T) 
a 2 L r 2 

' r 3 



2vk 2 ± 2vY' 



a -(l± a -\ 
vT \ k) 



(63) 



After expanding the expression (|63|) in power series in D 
up to linear terms we derive the main result (|1 7 |1 - (flfj f) . 



References 

1. P. Reimann, Phys. Rep. 361, (2002) 57-265. 

2. V. N. Smelyanskiy, M. I. Dykman, and B. Golding, Phys. 
Rev. Lett. 82, (1999) 3193-3197. 

3. J. Lehmann, P. Reimann, and P. Hanggi, Phys. Rev. Lett. 
84, (2000) 1639-1642. 

4. R. S. Maier and D. L. Stein, Phys. Rev. Lett. 86, (2001) 
3942-3945. 

5. J. E. Hirsch, B. A. Huberman, D. J. Scalapino, Phys. Rev. 
A 25, (1982) 519-532. 

6. I. Dayan, M. Gitterman, G. H. Weiss, Phys. Rev. A 46, 
(1992) 757-761. 

7. R. N. Mantegna and B. Spagnolo, Phys. Rev. Lett. 76, 
(1996) 563-566. 

8. R. N. Mantegna and B. Spagnolo, Int. J. Bifurcation and 
Chaos 8, (1998) 783-790. 

9. R. N. Mantegna and B. Spagnolo, Stochastic Processes in 
Physics, Chemistry and Biology, Lecture Notes in Physics 
(Eds. J. A. Freund and T. Poeschel, Springer- Verlag, 
Berlin 2000) 327-337. 

10. N. V. Agudov and B. Spagnolo, Phys. Rev. E 64, (2001) 
035102(R)-1 - 035102(R)-4. 

11. N. V. Agudov and B. Spagnolo, Stochastic and Chaotic Dy- 
namics in the Lakes (Eds. Broomhead D. S., et al, Amer- 
ican Institute of Physics, 2000) 272-277. 

12. N. V. Agudov, A. A. Dubkov, and B. Spagnolo, Physica A 
325, (2003) 144-151. 

13. A. Fiasconaro, D. Valenti, and B. Spagnolo, Physica A 
325, (2003) 136-143. 

14. A. N. Malakhov, A. L. Pankratov, Physica C 269, (1996) 
46-54. 

15. B. Spagnolo, A. A. Dubkov, and N. V. Agudov, Acta Phys- 
ica Polonica B 35, (2004) 1419-1436. 

16. R. Wackerbauer, Phys. Rev. E 58, (1998) 3036-3044. 

17. R. Wackerbauer, Phys. Rev. E 59, (1999) 2872-2879. 

18. A. Mielke, Phys. Rev. Lett. 84, (2000) 818-821. 

19. D. Dan, M. C. Mahato, and A. M. Jayannavar, Phys. Rev. 
E 60, (1999) 6421-6428. 

20. M. C. Mahato and A. M. Jayannavar, Mod. Phys. Lett. B 
11, (1997) 815-820. 

21. M. C. Mahato and A. M. Jayannavar, Physica A 248, 
(1998) 138-154. 

22. K. Matsumoto and J. Tsuda, Phys. Rev. A 32, (1985) 
1934-1937. 

23. M. Yoshimoto, Phys. Lett. A 312, (2003) 59-64. 

24. Xie Chong-Wei and Mei Dong-Cheng, Chin. Phys. Lett. 
20, (2003) 813-816. 

25. F. Apostolico, L. Gammaitoni, F. Marchesoni, and S. San- 
tucci, Phys. Rev. E 55, (1997) 36-39. 

26. C. R. Doering and J. C. Gadoua, Phys. Rev. Lett. 69, 
(1992) 2318-2321. 

27. M. Marchi, F. Marchesoni, L. Gammaitoni, E. Menichella- 
Saetta, and S. Santucci, Phys. Rev. E 54, (1996) 3479- 
3487. 

28. R. N. Mantegna and B. Spagnolo, Phys. Rev. Lett. 84, 
(2000) 3025-3028. 

29. P. Hanggi and P. Talkner, Phys. Rev. A 32, (1985) R1934- 
R1937. 

30. J. Masoliver, K. Lindenberg, and B. J. West, Phys. Rev. 
A 33, (1986) 2177-2180. 

31. M. A. Rodrigues and L. Pesquera, Phys. Rev. A 34, (1986) 
4532-4534. 



B. Spagnolo, A. A. Dubkov, N. V. Agudov: Enhancement of stability in randomly switching metastable state 9 



32. C. R. Doering, Phys. Rev. A 35, (1987) 3166-3167. 

33. V. Balakrishnan, C. Van den Broeck, and P. Hanggi, Phys. 
Rev. A 38, (1988) 4213-4222. 

34. U. Ziircher and C. R. Doering, Phys. Rev. E 47, (1993) 
3862-3869. 

35. J. Luczka, M. Niemiec, and E. Piotrowski, Phys. Lett. A 
167, (1992) 475-478. 

36. A. A. Dubkov, N. V. Agudov, and B. Spagnolo, Phys. Rev. 
E 69, (2004) 061103-1 - 061103-7. 

37. K. Furutsu, J. Res. Nat. Bur. Stand. 67D, (1963) 303-323. 

38. E. A. Novikov, Sov. Phys. JETP 20, (1965) 1290-1294. 

39. V. E. Shapiro and V. M. Loginov, Physica A 91, (1978) 
563-574. 



